HTL Resummation of the Thermodynamic Potential 



Andre Peshier 

Department of Physics, Brookhaven National Laboratory, Upton, 

NY 11973-5000, USA 



Abstract 

Starting from the ^-derivable approximation scheme at leading-loop order, the 
thermodynamic potential in a hot scalar theory, as well as in QED and QCD, 
is expressed in terms of hard thermal loop propagators. This nonperturbative 
approach is consistent with the leading-order perturbative results, ultraviolet 
finite, and, for gauge theories, explicitly gauge- invariant. For hot QCD it is 
argued that the resummed approximation is applicable in the large-coupling 
regime, down to almost twice the transition temperature. 
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1 Introduction 

One of the central issues of the ongoing heavy-ion program is the investigation of 
highly excited strongly interacting matter. As predicted by Quantum Chromody- 
namics (QCD), at energy densities of the order of IGeV/fm^ hadron matter will 
undergo a transition to a state of deconfined quarks and gluons. Despite the asymp- 
totic freedom of QCD, this quark-gluon plasma (QGP) is characterized by a large 
coupling strength in the regimes of physical interest. Therefore, nonperturbative 
approaches are required to describe this many-particle system reliably. 

This expectation has been demonstrated in the calculation of the thermodynamic 
potential of the hot QGP. The perturbative expansion, which is known through order 
0{g^) in the coupling shows no sign of convergence even orders of magnitude 
above the transition temperature ~ 170 MeV. Instead, with increasing order of the 
approximation, the results fluctuate more strongly as a function of the coupling (or 
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the temperature) along with a growing residual dependence on the renormalization 
scale. These features are not specific to QCD, but are also observed, e. g., in a scalar 
theory [Q, and are presumably related to the asymptotic nature of perturbative 
expansions. 

On the other hand, extensive studies in lattice QCD indicate for both the pure 
gauge plasma ||^, and for systems containing dynamical fermions ||^, that the ther- 
modynamic potential Q, scaled by the interaction-free limit, is a smoothly increasing 
function of the temperature, and approaches values ^ 80% of the ideal gas limit at 
T ^ 2Tc. This saturation-like behavior, as well as the rapid change of fl between Tc 
and 2Tc, has been interpreted within a quasiparticle picture P, assuming that for 
large coupling the relevant excitations of the plasma can be described by quasipar- 
ticles with effective masses depending on the coupling, as known in the perturbative 
regime. This phenomenological description amounts to a partial resummation of 
relevant contributions beyond the leading perturbative corrections. Therefore, the 
quantitative agreement of these models with the available lattice data is an indica- 
tion that resummations relying on an appropriate quasiparticle structure may indeed 
lead to improved approximations for the strong-coupling regime. 

For a scalar theory, this conjecture is supported by the so-called screened pertur- 
bation theory 0, where the conventional perturbative expansion is rearranged by 
adding and subtracting a mass term to the lagrangian and expanding in the massive 
propagator, treating the subtraction as an additional interaction. By relating the 
mass to the self-energy of the particles, important effects of the interaction are taken 
into account already at the leading order of the reorganized expansion, which indeed 
shows an improved behavior for large coupling ^ . The same idea of appropriately 
rearranging the lagrangian is applied in the hard thermal loop (HTL) perturbation 
theory Given the noteworthy properties of the HTL Green's functions which, 
in particular, satisfy fundamental sum rules and Ward identities in gauge theories, 
they arise as a preferable basis for reorganized expansions. In the context of ther- 
modynamics, the leading (zeroth) order contributions to the pressure have been 
calculated within the HTL perturbation theory for the hot QGP as well as for 
the degenerate plasma |jlO[- These nonperturbative expressions reproduce correctly 
the zeroth order of the conventional perturbation theory (i.e., the free limit). In 
addition, they already include effects of Landau damping and of electrical screen- 
ing, and hence reproduce at finite temperature the so-called plasmon effect of the 
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order 0{g^). However, they do not account correctly for the perturbative 0{g'^) 
contributions, since also the next-to-leading terms in the HTL perturbation the- 
ory contribute to that order. Moreover, as a formal point at issue, an ambiguous 
regularization scheme dependence related to uncompensated ultraviolet divergences 
which implicitly depend on the temperature or the chemical potential, will only be 
improved by the next-to-leading order calculation. This is similar to the situation 
for screened perturbation theory in the scalar case. 

A conceptually different approach is the selfconsistent ^-derivable approximation 
scheme |TT], which has been applied in for the scalar theory. At leading- loop 
order, this approximation is equivalent to the large- limit of the scalar 0{N)- 
symmetric model |]l3l, and leads to similar results as the screened perturbation 
theory, without having to drop the ambiguous temperature dependent divergences 
in the thermodynamic potential. In ||14| it was shown in QED that the entropy de- 
rived from the leading-loop ^-derivable thermodynamic potential can be expressed 
as a simple functional of the dressed propagators, which themselves have to be 
determined selfconsistently. This was generalized to QCD in Ref. [|15|, where a 
nonperturbative result for the QCD entropy has been obtained by approximately 
evaluating the entropy functional, at the expense of exact selfconsistency, with the 
HTL propagators. This approach avoids the nontrivial issues related to gauge in- 
variance and renormalization of the resummed propagators, and leads to a physical 
and formally well-defined approximation: it is gauge invariant, ultraviolet finite and 
reproduces the perturbative 0{g'^) result (the next-to- leading order is discussed in 
as well, see also the remarks below). Moreover, complemented with the two- loop 
running of the coupling strength, the HTL-resummed entropy matches the lattice 
results in the saturation-like regime, starting at temperatures T ^ 2Tc. 

From the entropy, given as a function of the temperature, the thermodynamic 
potential can be obtained, up to an integration constant. However, it is interesting 
to consider the thermodynamic potential itself in the ^-derivable approach. As a 
matter of principle, since in the framework of the HTL approximation (applied for 
the reasons mentioned above) the thermodynamical selfconsistency holds only ap- 
proximately, different calculations of the same quantity may lead to different results 
- which should be compared. More importantly, the approximate thermodynamic 
potential, expressed in terms of the self-energies which have to be determined as 
a function of the temperature, contains relevant information not given by the cor- 
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responding expression for the entropy. As will be shown for the scalar theory and 
conjectured for gauge theories, from the thermodynamic potential one can infer the 
range of validity of the approximation in the large-coupling regime. 

This paper is organized as follows. In section 2, the concept of ^-derivable 
approximations is resumed for the scalar theory. This provides, at leading-loop 
order, a solvable, yet representative, case of reference. Following |T^, selfconsistent 
and approximately selfconsistent results are derived and discussed with emphasis on 
the extrapolation to large coupling. For gauge theories, the approach is presented 
in section 3 for the Abelian case before considering QCD in section 4; parallels to 
the scalar theory are pointed out. The conclusions are summarized in section 5. 
Explicit expressions of relevant sum-integrals are relegated to the appendix. 



2 Scalar field theory 

2.1 ^-derivable approximations 

As an exact relation, the thermodynamic potential Q can be expressed in terms of 
the full propagator A by the (generalized) Luttinger-Ward representation |T6| , 

fi = ^tr[ln(-A-i)+An]-<l>[A], (1) 

where the trace is taken over all states of the system. The exact self-energy 11, which 
is related by Dyson's equation to the free propagator Aq and the full propagator by 

A-' = A,'-U, (2) 

can be represented diagrammatically as the series of the dressed one-particle irre- 
ducible graphs tin of order n. In the expression (|I]), the interaction-free contribution 
Qq is contained in the trace part, as is obvious from expanding ln(— A^^) + AH = 
ln(— Aq ^) + ln(l — AqII) + Aon/(l — AqII) in powers of 11, while the leading-order 
perturbative correction is entirely due to the $ contribution. The functional $ in 
(|l|) is given by the skeleton diagram expansion 

* = E^tr[An„]. (3) 

Thus, taking into consideration the combinatorial factors related to the number of 
propagators in each graph of this expansion, the self-energy is obtained diagram- 
matically from $ by opening one of the propagator lines in the individual graphs. 



4 



I.e., 

Consequently, the thermodynamic potential (P considered as a functional of A is 
stationary at the exact propagator defined by the solution of the implicit functional 
equation (H), 

sn 

This fundamental relation expresses the thermodynamical selfconsistency between 
microscopic and macroscopic properties of the system [ p!l[| . 

Commencing from this exact framework, selfconsistent approximations of the 



thermodynamic potential can be derived ITJ]: truncating the complete skeleton 
expansion of the functional $ at a certain loop order, and calculating selfconsis- 
tently the approximate self-energy analogously to (1), yields an approximation of 
the thermodynamic potential which is still stationary with respect to variations of 
the resummed propagator. In this ^-derivable approximation scheme, appropriate 
sets of diagrams contributing to the self-energy and to the thermodynamic potential 
are resummed, in such a way that thermodynamical properties can be calculated 
from the thermodynamic potential, utilizing thermodjTiamical relations, or directly 
from the Green's functions, which leads to the same, approximate result. 

2.2 Selfconsistent leading-loop resummation 

With the interaction of the scalar particles described by Ci = (— (7o/4!)0^, the 
functional $ and the corresponding self-energy are given at the leading-loop order 
in the skeleton expansion, with explicit symmetry factors, by 

= 3 ^ ^ 

n. = 12 =i2(^^)/(n,,T). (6) 

This truncation of the skeleton expansions is equivalent to considering the leading 
contributions of the expansion in the scalar 0(A^)-symmetric theory |T^, hence 



providing exact results in the limit — ^ oo. In terms of the conventional pertur- 
bation theory, it amounts to a complete resummation of the so-called super-daisy 
diagrams [0. 
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Since in the 2-point contribution is local, 11;; is just a mass term, which 
reduces the implicit functional equation to a transcendental, yet nontrivial, relation. 
In the imaginary time formalism, regularizing the spatial momentum integrals in 
(i = 3 — 2e dimensions in the MS scheme, the trace over the momentum K = [ko, k) 
is defined as 

with the renormalization scale /2 and Euler's constant 7, and where ko = i 2mTT are 
the bosonic Matsubara frequencies. The function 

/(M^ T) = ^ = I'{M') + I\M\ T) (7) 

is decomposed into two parts: a contribution P which is not explicitly temperature 
dependent (below, M will depend on T), and associated to the quantum fluctuations 
of the vacuum, and a part due to the thermal medium, which are both given in 
the appendix. While the thermal fluctuations are cut off by the Bose distribution 
function nb{x) = [exp(x) — the ultraviolet divergence of the vacuum contri- 
bution is apparent in a pole term oc M^/e. In the gap equation (||), identifying 
with the mass term rriQ + 11;; of the propagator A;;, this divergence is tempera- 
ture dependent. The contribution oc rriQ/e is absorbed in the physical mass by the 
conventional mass renormalizationQ. Focusing, to simplify the discussion, on the 
temperature dependent part oc Hu/e, only the case of massless particles (as relevant 
for the ultrarelativistic gauge plasmas studied below) is considered in the following, 
with Ao\ko, k)=K^ = kl- 

The thermal divergence in (^ requires the renormalization of the bare coupling 
Qq. The physical coupling can be related to the vacuum scattering amplitude. In 
the approximation considered, resumming the set of chain diagrams, 

X = X + i2p<. 

the renormalized coupling at the momentum scale s = is determined by 




L{P') = [ 

JK 



(P - Ky 167r2 



"'^In the present approach, the mass counter term is obtained from a 'counter loop' contribution 
(5$ — itr A((5m)^ to $, which does not affect the thermodynamic potential. 
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Omitted here are the crossed diagrams in the scattering amphtude; these induce 
graphs with a different topology in the expansion of the thermodynamic potential 
and the self-energyQ, 



Expressed in terms of the renormalized coupling, the Dyson equation reads 

(9) 



The right hand side of this implicit gap equation for is finite and independent 
of fi. Describing the system under consideration by the value of g'^{s) at a specific 
scale s, and taking into account that 9~\s') = g-\s) - In(s7s)/(327r2) according 
to eqn. the solution of @ is invariant under the rescaling s — > s' . In weak 
coupling, by expanding I'^iliu.T) in small II/z/T^, the resummed leading loop self- 
energy reproduces the perturbative result 



n -»»^ 



(10) 



up to next-to-leading order, where the coupling is still unrenormalized. In the gap 
equation (H), the selfconsistent resummation leads to a nontrivial interplay between 
vacuum and thermal fluctuations. In particular, as discussed in detail in [0, for 



not too large coupling the gap equation has two solutions, see Fig. |I|. While the 
smaller one is associated with the perturbative approximation, the second, tachyonic, 
solution is exponentially large for small g'^ and thus of no physical relevance. For the 
choice —s = and couplings larger than the value g^^^^ ~ 100, the gap equation 
has no solution. Clearly, the present approximation is physically justifiable only 
for couplings below g'^^^, where both solutions of the gap equation are not of the 
same order as the maximal value 11™^^ ~ 4.2 T^. The perturbative result ([lOD , on 



the other hand, shows the typical features of asymptotic expansions: with higher 



^ I. e., the contribution of the set of super-daisy diagrams to the self-energy is finite after renor- 
mahzing the couphng in the same class of graphs. As a consequence of the present approximation, 
therefore, the running of is determined by a /3 function which differs from the perturbative result 
by a factor of 1/3. It agrees with the f3 function of the scalar 0{N) model in the large- limit 
considered in [l^ , where the crossed diagrams in the scattering amplitude are suppressed. Hence, 

as defined in (0) is related to the = 1 running coupling only at a single value of the scale. 
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Figure 1: The self-energy, in the scalar theory, as a function of the coupling. 
Compared are the selfconsistent (sc) solution of the gap equation, with the 
coupling renormalized at the scale —s = T^, and the leading (l.o.) and next- 
to-leading order (n. l.o.) perturbative results depending on the bare coupling. 



order, the accuracy of the approximation is improved only in a limited range, which 
becomes smaller with increasing order, of the expansion parameter. 

Having solved the gap equation for a given value of g'^{s), the approximation Qa 
of the thermodynamic potential can be calculated as a function of 11;/. From the 
vacuum part J° of the function 

J{M\T) = ln(-Ao 1 + K'P) = + J^{M^, T) , (11) 

which is given explicitly in the appendix, it appears that the contribution 

(fi + <|.), = l^[ln(-A,i)+A,n, 



2- 

+ finite terms (12) 



e e 



647r2 

to the thermodynamic potential contains thermal divergences which are only partly 
compensated between the ln(— A^^) and the AH term, as indicated by the bracket. 
In fact, they are cancelled by the remaining $ contribution. With the two-loop $ 
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functional evaluated with the selfconsistent propagator, 

<^u[^a] = l^^iina, (13) 

where the renormalization of the coupling has been taken into account, the re- 
summed thermodynamic potential takes the form 



(14) 



which is, indeed, finite and independent of the renormalization scale. In terms of 
the functions (0) and (0), the pressure p = —fl (the volume is set to unity) is, in 
the selfconsistent approximation, 

pu = -J^iUu, T) - 1 n,, I^iUu, T) + ^ . (15) 

The first term is the pressure of a free gas of quasiparticles with mass n|/^, while 
the other terms represent the interactions among these quasiparticles. The last 
term is only implicitly temperature dependent. It stems from the vacuum part of 
the resummed contributions, and leads at larger 11/; to an increasing ratio of pa to 
the free pressure po = as shown in Fig. ^. This behavior, however, is 

physically not relevant. Instead, it indicates a breakdown of the approximation: as 
a result of the stationary condition (^ of the thermodynamic potential, the value of 
the self-energy at the point where pu has its minimum coincides with the maximal 
solution n™'^^ of the Dyson equation, 

nr^~4.2T2, pr=P«(nr^)~0.73po. (16) 



Therefore, the behavior of the thermodynamic potential ([1^) provides the same 
strong criterion for the applicability of the approximation as that of the solution of 
the gap equation, which, in both cases, is due to the interplay between the resummed 
contributions of the thermal and the vacuum fluctuations. 

It remains to note that similarly to the expression (^ for the dressed self-energy, 
the approximation ([T5|) of the pressure, which resums terms of all orders in the 
coupling, agrees with the perturbative result 



1 _ ^ + 15 



87r2 4! 27r3 




(17) 
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Figure 2: The selfconsistent pressure and the entropy (in units of the free 
values) in the scalar theory, as functions of the resummed self-energy. These 
results coincide with the HTL approximations of p and s as functions of 11*. 
Physical relevance can be attributed to these approximations only in the regime 
where p is a decreasing function of 11, see text. Also shown for comparison are 
the leading and next-to-leading order perturbative results (for both p and s) 
depending on Qq/AI = U* /T"^. 



up to next-to-leading order. As asymptotic expansions, the perturbative results 
start to fluctuate with increasing order, and the range of predictability is rather 
small already for the next-to-leading order expression, as apparent in Fig. ^. The 
resummed approximation, on the other hand, continues smoothly into the large- 
coupling regime, as expected on physical grounds for the exact result. 

Shown also in Fig. H is the entropy (density) related to the pressure by the 
thermodynamic relation s = —dQ/dT = dp/dT. In the leading-loop approximation, 
from eqn. (0), 



sii = - 



dT 



4 



dT 



dpii 



dUg 
dT 



the second and the third terms cancel for the selfconsistent solution of the gap 



10 



equation, hence 

dr{IiiuT) 



su 



dT 



{li 

Till 



Contrary to the pressure (H), the entropy su, as a measure of the population of the 

1/2 

phase space, is equivalent to that of a system of free quasiparticles with mass 11;/ 
|l^ . Therefore, su is a monotonically decreasing function of the self-energy, without 
any indication of the breakdown of the approximation. However, according to the 
preceding considerations, only values of the entropy larger than s™™ ~ 0.60 sq are 
justified in the present approach. 

Alternatively (not specific to the scalar theory), the entropy derived in the self- 
consistent approach with the two-loop approximation of $ can also be evaluated 
from a functional, independent of $, of the dressed propagator |T^. This direct 
calculation benefits from the observation that the entropy is, in contrast to the pres- 
sure, sensitive only to the thermal excitations of the system and, thus, a manifestly 
ultraviolet finite quantity. Formally, this is related to the fact that the function 
piiillii) indicates a breakdown of the approximation, but suijlu) does not. This 
statement does not contradict the thermodynamical relation between the pressure 
and the entropy since, to reconstruct the pressure pu from eqn. ([T8|) , the selfconsis- 
tent self-energy has to be known as a function of the temperature. In that sense, the 
pressure expressed in terms of the self-energy provides relevant information not con- 
tained in the entropy - which makes the pressure an interesting quantity to consider 
for cases where the self-energy cannot be resummed self consistently, as examined in 
the following. 



2.3 Approximately selfconsistent resummations 

With regards to gauge theories considered below, it is instructive to study an addi- 
tional approximation within the leading-loop <l>-derivable approach, by solving the 
gap equation only approximately. 

The leading-order perturbative contribution of the resummed self-energy is 
obtained from the tadpole graph with the bare propagator and agrees with the HTL 
approximation (denoted by a star), 

= ^ . (19) 
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At this level of approximation, no renormalization is required (the perturbative 
vacuum divergence vanishes in dimensional regularization) . As put forward in ||I5| , 



a resummed approximation of the entropy can be obtained from the aforementioned 
entropy functional evaluated with the dressed propagator 

A* = [Ao' -n*]"\ (20) 

making use of the observation that this functional does not dependent on $;/ for 
any ansatz of the propagator. The resulting approximation of the entropy is then 
given again by the free quasiparticle expression, and the only difference to the self- 
consistently resummed approximation ([l8|) is the quasiparticle mass, which is now 
determined by (|19|). In contrast to the selfconsistent result, however, this approxi- 
mation reproduces the perturbative series only to order O((7o)- The next-to-leading 
correction is underestimated by a factor of 1/4 since, diagrammatically, the set of 
daisy (ring) diagrams is included correctly only after the resummation of the self- 
energy. Indeed, replacing 11* by the next-to-leading order perturbative expression 
(|TUp, as considered in |TB[, yields the correct O(5'o) term for the entropy. It is 
recalled, however, that the 0{gQ) contribution to the self-energy has its origin in 
the screening of the dressed propagator, and that the related thermal divergence 
~ g^^goTy, which cancels in the selfconsistent gap equation, has to be dropped 
here in the spirit of perturbation theory^. 

Returning now to the thermodynamic potential, one could approximately evalu- 
ate the functional (|1]) with the two-loop contribution to $ and the HTL propagator 
(pOl). This approximation indeed resums the set of the individually infrared diver- 
gent daisy diagrams, and thus reproduces the perturbative result up to the order 
0{gQ) for the plasmon term. However, $;;[A*] does not combine with the second 
term of | tr[ln(— A* ^^) -|- A*n*] as for the selfconsistent approximation, so the total 
expression contains a temperature dependent divergence oc (11*)^ oc (goT)^, and is 
therefore relevant only as a perturbative expansion up to order 0{gQ). 

Alternatively, one can approximate the $ contribution by 



l^A*n*, (21) 



■^On the other hand, to obtain improved quantitative results for the entropy at larger coupling 
strength, the authors of Ref. used a Fade approximation of the next-to-leading order self- 
energy, which contains higher powers of go. 
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which resembles the selfconsistent contribution in (|T3|). As required for an appro- 
priate approximation (since the leading-order contribution to the thermodynamic 
potential arises from the $ functional), $* agrees to order OIqq) with ^u. Further- 
more, the resulting approximation for the thermodynamic potential. 



:22) 



is obviously ultraviolet finite because it has the same functional dependence on 11* 
as the selfconsistent approximation The expansion of the pressure for small 

U*/T^ = g^/Al is 

P^ = Po [l-0]-[2-l] — - + [1-3/4] — (-) +... , (23) 



(the brackets indicate the contributions of the two terms in eqn. (|22|) ). This re- 
produces the perturbative result (|T^ to order C(fi'o)) but underestimates the O^g^) 
term by a factor of 1/4 , for the same reason and with the same implications as given 
above for the entropy. Even so, the approximation (^2]) is physically significant, since 
the self-energy (in general, the mass scale of the self-energy) is in principle a mea- 
surable quantity. It is more important, however, that it also provides a restricting 
criterion for the validity of the approach to approximate the thermodynamic po- 
tential in terms of the HTL propagator. Clearly, the approximation f2*(n*) cannot 
be justified in a regime of strong coupling where even the selfconsistent approxima- 
tion is not physical. In other words, the minimum of p*{Il*) provides, even without 
having at hand the selfconsistent solution of the gap equation, a strong limit of ap- 
plicability of the approximately selfconsistent approach. This fact will be important 
in the following discussions of gauge theories. 



3 QED 

In this section, an ultrarelativistic electron-positron plasma, described by QED, is 
considered at temperatures much larger than the electron mass. 

For gauge theories, an additional requirement for a formally consistent approxi- 
mation of a physical quantity is gauge invariance. In the ^-derivable approach (un- 
less solved exactly), by resumming thermodynamically selfconsistent sets of graphs, 
this requirement is, in general, not satisfied, since the two-point functions are dis- 
tinguished in the hierarchy of Green's functions. Besides, it is not obvious how a 
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renormalization of the coupling constant can be accomplished technically, to account 
for the thermal divergences in the resummed Dyson equations. However, it will turn 
out that the leading contributions to the resummed thermodynamic potential arise 
as in the scalar theory from the HTL parts of the propagators, which are, in fact, 
gauge invariant, and renormalized by the usual vacuum counter-terms. Therefore, 
an explicitly gauge-independent nonperturbative approximation of fl can be derived, 
which exhibits analogous features as the approximately selfconsistent resummation 
in the scalar case. Later I will argue that this approximation allows one to conjecture 
about the large-coupling behavior of the leading-loop resummation. 

The exact thermodynamic potential can be expressed as a functional of the 
photon propagator D and the electron propagator S, which are by Dyson's equation 
related to the respective self-energies 11 and S, a^ 

n = ^ Tr[ln(-D-^) + DU] - Ti[\n{-S-^) + SE] - S] . (24) 

In the boson contribution, the trace is taken over the four-momentum as well as 
over the Lorentz structure, while in the fermion part the trace includes the spinor 
indices. The functional $[-D, S] is given by the series of skeleton graphs with exact 
propagators. It is related to the self-energies by 

(5$ (5<I> 

so the thermodynamic potential is stationary at the full propagators. 

By the projectors Pj, = g^, - K^KJK^ - P^, and Vj:, = -K^KJK^ where 
K = [K(Ku) —uK'^]/[{Kuy — K'^Y^'^ and u is the medium four- velocity, the inverse 
photon propagator is decomposed into the transverse (T) and the longitudinal (L) 
part as well as the covariant gauge-fixing term, 

i=T,L ^ 

Introducing the 'projectors' V± = \ {IjC^ljC) (the index denotes the ratio of chirality 
to helicity), the fermion propagator can be written in a similar way as 

i=± 



''Depending on the gauge, the contribution of the AbeUan ghost fields, which compensate the 
non-transverse degrees of freedom, is included implicitly. 
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In terms of the scalar propagators Aj, with the degeneracy factors dr 
and d± = {d + l)/2, the first two contributions in eqn. (p^ read 



d — 1, di = I 



^] + $ 



^ ^ I E^^^ [in(-Ari) + A,n,] - ln(-Ao-^) | 



(26) 



where the subtractive contribution of the Abehan ghost fields, which otherwise de- 
couple, is included in the boson term. 

To leading-loop order, ignoring for a moment the question of gauge dependence, 
the $ functional and the self-energies are determined by 



n, 



(27) 



In the following, as motivated above, the self-energies are approximated by their 
HTL contributions. This approximation, as a matter of fact, undermines the self- 
consistency of the ^-derivable approach. As discussed for the scalar theory, without 
the resummation of the self-energies in (|27|) , the set of ring diagrams contributing to 
the thermodynamic potential at next-to-leading order is included only incompletely, 
and the coupling constant e remains bare. Even more important, it is not obvious a 
priori that a physically meaningful approximation of Q can be formulated in terms 
of the HTL propagators. These are derived for soft momenta, much smaller than 
the temperature, whereas thermodynamics is sensitive to the momentum scale T. 
It will turn out, however, that the dominant contributions originate from the vicin- 
ity of the light cone, where the HTL approximation is appropriate, even at large 
momenta. 

In the HTL approximation, the self-energies are given by ^ 



- m; ± s 

2 " 



-2n, fl{ko,k) 



ko ko-k 
1 H r In 

2k ko + k 



M2 

S(A;o,fc) = ^ 



ko K'^ ko-k 
m 

k 2k^ ko + k 



The quantities 



Mi 



Mt 



i2t 



(29) 
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are referred to as the asymptotic masses (squared) of the transverse photon and the 
electron particle excitations, respectively, since their dispersion relations approach 
mass shells for momenta k ^ eT. The longitudinal photon (plasmon) mode and 
the electron hole (plasmino) excitation, on the other hand, possess an exponentially 
vanishing spectral strength for k ^ eT when approaching the light cone. 

The contribution ( P^D to the thermodynamic potential can be approximated by 
evaluating the functional with the HTL propagators (see appendix and, for details 
of the calculation, p], [l^). It is worth pointing out that the ln(— A~^) and the 
AjAg ^ contributions of the fermions, summed over i = ±, are ultraviolet finite in- 
dividually, while in the respective boson contributions only the most severe thermal 
divergences oc M^/e^ cancel between the transverse and longitudinal terms. Sur- 
prisingly, although now the self-energies are nonlocal and have imaginary parts, the 
structure of the overall divergence is in an apparent analogy to the scalar case (pr2D, 

where the terms in the bracket originate from the first and second contribution in 
iTr[ln(— D^^) +Dn]. Again, this divergence is temperature dependent and, thus, 
expected to be cancelled by the remaining $ contribution. 

In a selfconsistent approach, the functional evaluated with the selfconsistent 
propagators could be expressed as traces over the self-energies, 

*„ = l™ = -iT.S„E,„ (31) 



1 2 
- + - 

e e 



+ finite terms , (30) 



analogously to the identity ( P^ ) in the scalar theory. However, the naive replacement 
of selfconsistent propagators and self-energies by their HTL approximations, as in 
(pTD for the scalar case, leads to different results for the bosonic and the fermionic 
traces even at order 0{e^), since the expressions are dominated by hard momenta 
where contributions neglected in the HTL approximation become important. Even 
so, there is an unique combination of the two traces which keeps track of these terms. 
Denoting the photon momentum in | Q by and the fermion momenta by Qi,2, 
this diagram (with bare propagators for the 0{e^) contribution) can be represented 
as a double sum- integral over an expression with a numerator = ^ Q\ — Q\- 
Closing the external legs of the boson self-energy in the HTL approximation amounts 
to neglecting the term in A^. Tracing over the fermion HTL self-energy, on the 
other hand, neglects one of the terms. Thus, all terms are accounted for twice 
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in the sum over all three possibilities to approximate one of the momenta as soft. 
Accordingly, the specific combination 

= ^ TiD*n* - - Tr5*S* (32) 
4 2 ^ ^ 

is equivalent, at leading order, to the $ contribution (|31|). As shown in the following, 
this unique approximation indeed leads to a well-defined resummed approximation 
of the thermodynamic potential. 

It is first emphasized that the complete expression resulting from the HTL ap- 
proximation of eqn. ( p6|) and (|3^) , which can be written in a compact form as 



= -Tt 
2 



ln(-D*-i) + ^D*U* 



Tr 



ln(-5'^ -1) + i 



(33) 



is analogous to the corresponding expression (^2]) in the scalar theory. After the 
precedent discussion of the structure of the divergences of the individual terms in 
the contribution (pOD, this implies that f2* is an ultraviolet finite quantity and explic- 
itly independent of the regularization scale. In the interaction- free limit, it reduces 
to the thermodynamic potential ^Iq = — ((ir + 1 2(i±) |^ of an ideal gas of photons, 
with the non-transverse degrees of freedom compensated by the ghost contribution, 
and electrons/positrons. Moreover, the leading-order perturbative result is repro- 
duced since it originates, as noted above, solely from the $ contribution. From the 
individual terms in (|33|), the boson and fermion contributions of order O(e^) arise 
in analogy to ( P5| ) in the scalar theory, 

= ^ ([2 - l]M'^ + [2 - l]Mf ) = fifr . (34) 

As in the calculation of the HTL entropy, the leading-order contribution is 
therefore determined entirely by the behavior of the self-energies near the light cone, 
where the HTL approximation is, in fact, valid also asymptotically [T^, which justi- 
fies a posteriori the usage of this approximation in the present approach^. However, 
in contrast to the entropy which is manifestly ultraviolet finite, the agreement of 
Q* with the perturbative result is directly related to the cancellation of the thermal 
divergences. This aspect provides another, formal, argument for the approxima- 
tion (|3^) of the $ contribution: any other linear combination of traced boson and 



^This also supports the phenomenological quasiparticle models applied (for QCD) in which 
consider the thermodynamically relevant excitations as a gas of quasiparticles with effective masses 
equivalent to the asymptotic masses. 
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fermion self-energy contributions would result in either an uncompensated thermal 
divergence or an incorrect perturbative limit. 

The next-to-leading order term of the perturbative expansion, as already dis- 
cussed, cannot be expected to be reproduced in the present framework. Neverthe- 
less, it is a noteworthy remark that 

^^:lo = -[l-3/4]^(2M^T/^ (35) 

which originates from the static longitudinal parts of the \n{~D~^) and the DU 
contribution in (p3|), as indicated by the bracket, underestimates the perturbative 
result f^nfo* = — T/(127r)(e^T^/3)'^/^ again by a factor of 1/4, for the same combi- 
natorial reason as in the scalar theory, cf. (pSf ). As for the scalar model, the 0{e^) 
term could be reproduced correctly by evaluating the two-loop $ diagram with the 
dressed HTL propagators. Not surprisingly, this kind of approximation results also 
in the present case in an uncompensated thermal divergence oc (eT)^ and is, thus, 
relevant only as an expansion up to next-to-leading order. 



The full expression (33) resums terms of all orders in the coupling constant and 
has to be evaluated numerically. The pressure p* = —Q* turns out to be a mono- 
tonically decreasing function of the coupling |T^, with the anticipated behavior of a 
physical approximation: it yields a smooth extrapolation of the leading-order pertur- 
bative result to the large- coupling regime, where it is enclosed between the leading 
and the next-to-leading order approximations which, as asymptotic expansions, pre- 
sumably represent lower and upper boundaries for the exact result. For large values 
of e, the behavior of the pressure p* is determined by terms oc M^^ oc (eT)^ which 
stem, as shown in the scalar theory, from the interplay between the quantum and 
thermal fluctuation^. The fermion contribution oc is negative and overcompen- 
sates the positive photon contribution, and the total pressure p* becomes negative 
when extrapolated beyond e ~ 8.5. While in this regime the resummed approxi- 
mation is certainly not applicable, it will be argued in the following section that 
the boson contribution indicates the breakdown of the approximation already for a 
smaller coupling strength, at a = e^/(47r) ~ 1.5. 



In the present case, parts of these terms are imphcit in the individual quasiparticlc and Landau- 
damping contributions, depending on the chosen subtraction terms, see Appendix. 
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4 QCD 



In this section, I consider QCD, the original problem of interest. Within the frame- 
work of the HTL approximation, the quark propagator is analogous to the electron 
propagator in the Abelian plasma. Accordingly, the contribution of the quarks to 
the thermodynamic potential is analogous to the electron contribution in QED, so 
the following considerations focus on the gauge sector of QCD. 

For a pure gauge plasma, the thermodynamic potential is expressed as a func- 
tional of the gluon propagator D = (-Dq ^ — H)^^ and, in covariant gauges, the 
propagator G = {Gq^ — S)^^ of the anticommuting bosonic ghost fields, by 
1 



n 



Tr[ln(-D-^) + DU] - Tr[ln(-G"^) + GH] - $ 



(36) 



where the summation over the color indices is implicit in the trace. With the same 
remark about gauge invariance as in the previous section, the functional $ and the 
related self-energies are given to leading-loop order by 



n. 




H — araim 
2 \ / 



-|- OTi 



(37) 

In contrast to the photon polarization function, the gluon self-energy need not be 
transverse and it depends, in general, on four scalar functions. However, approx- 
imating the selfconsistent solutions of the coupled Dyson equations by their HTL 
contributions, as in QED, the gluon self-energy coincides with the Abelian expression 
(pSD up to the replacement of by the asymptotic gluon mass 



Mi 



2rr2 



6 



(38) 



with iVc = 3 for QCD. The ghost self-energy vanishes in the HTL approximation. 
Therefore, after rescaling the asymptotic mass and taking into account the color 
degrees of freedom, the contribution 

1 



{VL + $)* = - Tr[ln(-D* -^) + D*H*] - Tr[ln(-G^^)] 



(39) 
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to the thermodynamic potential is equivalent to the respective photon-ghost contri- 
bution in QED. The $ functional, however, has a more complicated topology than 
its Abelian counterpart, and requires a more detailed discussion. 

Evaluated with the selfconsistent solutions of the coupled Dyson equations, the 
contributions of the individual graphs for could be expressed as 



4 



"2 TrGji^i/ 



(40) 



where, in obvious notation, the terms 11^; denote the individual contributions to 11/; 
in (pTl). Note that at this level of approximation, unlike in QED, $// cannot be 
represented as a linear combination of the gluon and ghost self-energies traced over 
their external momenta. 

As for QED, the $ contribution cannot be approximated by naively evaluating 
the expressions (|40|) within the HTL approximation, since the terms of higher order 
in the external momentum, which are neglected in the HTL approximation of the 
self-energies, are important in the traces over the whole phase space. Again, to keep 
track of these terms, the leading-order contributions of the individual expressions in 
(PUD, which are obtained by replacing the dressed propagators by bare ones, are now 
considered explicitly. The first term can be represented as a double sum- integral 
over an expression with a numerator = + |-P^. Closing the external legs of 
the HTL contribution to H^s(P) amounts to neglecting the term ^P^ in A^, thus 

6 " " lo 2 6 lo 
The leading-order contribution of the diagram containing the gluon tadpole is sat- 
urated by the HTL contribution. 



(41) 



1 



lo 



(42) 



lo 4 

Finally, it is noted that the contributions omitted in the HTL approximation of the 
last two terms in (|40|) are taken into account in their sum, hence 



1 
4 



- TrZ^«Hf 



lo 



- TrG//^/; 



lo 



4 

- TrD*H§^* 
4 



lo 



lo 



- TrG^H' 
4 



lo 



(43) 
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since the ghost self-energy vanishes in the HTL approximation. Combining the right 
hand sides of the expressions (PD-(^), the $ contribution is approximated by 



= ^ Tr[D*(n3g* + n^g* + W^*)] = ^ TtD*U'' , (44) 



which is gauge independent. As for the corresponding expression (|32D in QED, the 
compact form of (0) is not a surprise. Rather, it reflects the fact that the combi- 
natorial factors related to the number of propagator lines in the skeleton diagrams 
also account for how many of them can be soft in the HTL approximation. 

The resummed approximation of the SU (Nc) thermodynamic potential reads 



Q* = -Tr 
2 



ln(-Z}^-i) + ^D*n^ 



Tr[ln(-Go^)] (45) 
which, as already stated for the contribution (|39D, resembles the photon-ghost part 



in QED in eqn. (|33D . Hence, the analogy of the photon and the gluon propagators 
in the HTL approximation becomes also apparent in the thermodynamical proper- 
ties, as already anticipated for the fermion contribution. Therefore, the properties 
emphasized for the approximation (|3^) in QED hold also for (^5]) in QCD, namely, 
the parallels to the scalar case and the structure of the thermal divergences (which 
arise solely from the boson contribution). Accordingly, the expansion of (|45| ) in 
Mg/T^ = I A^cS'^ reproduces the leading-order term of the perturbative result 



Qpert ^ _2/jy2 _ 

^ " ^90 



\ 15 N^g' 15 f^^N^g'V^" 
87r2 6 47r3 I 6 / 



(46) 



but underestimates the next-to- leading correction by a factor of 1/4, as in the scalar 
theoryQ. The resummed pressure p* = —Q*, evaluated numerically from eqn. (|45| ) 
and depicted in Fig. ^ as a function of the mass scale of the gluon self-energy, dis- 
plays a striking resemblance to the scalar case shown in Fig. |^. At small coupling it 
matches the leading-order perturbative result, but it decreases less fast with increas- 
ing Mg. For larger coupling, the pressure is dominated by a contribution oc M^, 
which is positive as in the scalar theory and leads, eventually, to an increasing 
behavior of the approximation, with a minimum at 

M|~3.1T2, p*^,^ = p\M^g)r^ 0.82 po. (47) 



''Replacing Mg Mg/[1 + ^{Nc/Sy^'^g], as suggested in reproduces the next-to-leading 
order term in 0^'^''* and might yield an improved resummed approximation for f2 at larger coupling 
strength, similar as in the scalar case. 
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Figure 3: The pressure p*, from eqn. (|45|) , and the entropy shtl [Hi of the 
pure gauge plasma in the HTL approximation (scaled by the free values) as 
functions of the asymptotic mass. Depicted as well are the leading and next- 
to-leading order perturbative results as functions of \Ncg'^ = Mg/T'^. 



These values compare to the numbers given in (|T6D for the scalar theory, which 
indicated in this simpler case, where the self-energy is just a mass term, that the 
approximation breaks down. Given the similarities of the expressions (^) and 
(PSD, the approximation for the thermodynamic potential in the SU{Nc) theory, 
with a much more complex structure of the propagators, cannot be expected to 
be appropriate for couplings corresponding to asymptotic gluon masses larger than 
Mg. The entropy shtl calculated in [0 from the leading loop entropy functional 
with the HTL propagators and depicted also in Fig. |^, on the other hand, shows no 
indication of a breakdown of the approximation at Mg. Drawing the parallel to the 
scalar theory, however, the extrapolation of shtl cannot be physically meaningful 
beyond the point where the approximation of the pressure breaks down. 

Taking an optimist's point of view, one may hope that the HTL-resummed pres- 
sure represents a reasonable approximation for 0.8 po ~ < Po, which, according to 
QCD lattice simulations & ffl, corresponds to temperatures larger than T ~ 2.5 Tc. 
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This conjecture is supported by the observation that the HTL entropy, supplemented 
('by hand') with the two-loop running coupling, is systematically off the SU{3) lat- 
tice results P below 2Tc, but starts to match the data for larger temperatures, 
just at T. Therefore, the HTL approximations are presumably predictive even for 
moderate large coupling, where the conventional perturbative results are no longer 
meaningful. This allows to systematically address physically relevant questions, e. g., 
related to finite chemical potential^, which can not yet be answered by QCD lattice 
calculations. 

For even larger coupling, in the vicinity of the confinement transition, where 
finite-temperature lattice simulations predict a rapid change of the thermodynamic 
potential, more sophisticated approximations are needed. These may also clarify 
the question whether the qualitative change in the behavior of the thermodynamic 
potential, at about 2.5 Tc, is related to the breakdown of the leading-loop approxi- 
mation at T, which in the present approach remains a speculation. 

5 Conclusions 

In this paper, nonperturbative resummations of the thermodynamic potential have 
been considered, starting from the ^-derivable approximation scheme at leading-loop 
order. By approximating the dressed propagators by their HTL contributions, at the 
expense of strict thermodynamical selfconsistency, physical and formally well-defined 
results have been obtained for the hot scalar theory, as well as for QED and QCD: 
they are ultraviolet finite, agree with the leading-order perturbative results and 
are, for gauge theories, explicitly gauge invariant, as are the approximations for the 
HTL-resummed entropy [0. I stress, however, that the finiteness the approximate 
pressure is, due to cancellations of ultraviolet divergences of individual terms, a 
more subtle issue which does, on the other hand, provide a formal argument for the 
present approach. 

Compared to the perturbative predictions, the HTL-resummed approximations 
display an improved behavior when extrapolated to larger coupling strength. For 

^By Maxwell's relation, as pointed out in |20| for the phcnomenological quasiparticle models 
the thermodynamical properties at finite temperature and at finite chemical potential are closely 
related, which allows us to infer, from the available finite-temperature lattice data, the equation 
of state of the QGP at nonzero baryon number. 
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all cases considered, at very large coupling strength, the resummed pressure is dom- 
inated by terms oc M^, where the mass scale M of the self-energies depends on 
the temperature and the coupling. These terms stem from the vacuum parts of the 
resummed contributions to the thermodynamic potential and lead, for bosons, to 
a pressure which increases as the coupling does. It has been shown for the scalar 
theory, and argued for gauge theories, that this feature is an indication of the break- 
down of the leading-loop resummed approximation. 

In QCD, this happens when the pressure is about 80% of the free limit, at a 
temperature T ~ 2.5 Tc. This explains the fact that the HTL-resummed entropy [ITH 
starts to deviate systematically from the lattice results below 2.5 Tg. On the other 
hand, the HTL approximation matches the lattice data for T ^ T. Therefore, the 
resummed approximations can be expected to be physical for rather large coupling, 
in a nonperturbative regime where the lattice simulations predict a saturation-like 
behavior of the thermodynamic potential. For even larger coupling, however, more 
sophisticated approximations are required for a detailed theoretical understanding of 
the thermodynamics of the QGP in the close vicinity of the confinement transition. 



A Appendix: The sum-integrals 



In this appendix, the relevant boson sum-integrals are given explicitly in the MS 
regularization scheme in = 3 — 2e spatial dimensions. 

Considering first the scalar theory, the vacuum and the thermal contributions of 
the trace /(M^, T) of the free propagator with mass M, defined in (|^), are 



J°(M2) 



16^ 



d?k ribiuJu/T) 



(2vr)= 



where ujm = {k"^ + M^)^/^. The contributions of the function J = \ tr ln(M^ — K"^) 
introduced in (0), which is related to / by differentiation with respect to M^, read 



J^{M^,T) 



e M2 



(49) 
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For QED and QCD, apart from the so-called quasiparticle contributions originating 
from the poles uji{k) of the propagators, there are Landau-damping contributions 
from the imaginary parts of the self-energies, i. e., from below the light cone in the 
HTL approximation. The individual contributions are regularized by appropriate 
subtraction terms. The expressions given below apply for both gauge theories (in 
QCD, the color trace yields just a factor), with Mf, being the respective asymptotic 
boson mass. The expression 



1 *N 



i=T,L,gh 

E di 

i=T,L,gh 



fc3 



+ Tln(l - exp[-cu,/T]) + / — (1 + 2n 

Jo /TT 



AT 



327r2 



I 13 TT^ 2 4 ,2 

n-^ + \n2 + -\-n^2 

e Mi 9 3 9 3 



(50) 



is analogous to the function J in the scalar theory, cf. (^91). In particular, the relevant 
number of degrees of freedom is apparent in both the thermal (quasiparticle) and 
in the divergent vacuum contribution. In (pUD, the contribution of the ghost fields 
is, formally, included as another degree of freedom with dg^ = —1 and Ugh{k) = k. 
The angles (pi are defined as (f)gh = and (f)T,L = Discln(— A^^/), and the infrared 
finite subtraction terms were chosen as Agh = and 



E d^At 



—di 



i=T,L 



Mi 



k 

2+4F + 



Ml 



/3_ A{k^ + Mi) - 

8k{k^ + Mi) [2 



Mi 



doj 
2^ 



Imn 
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1 + 



Ren 
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Refi A;2 ^ 

k^ + Mi 



The contribution resembling ^ M^J in the scalar theory is 

1 (1 + 2nb)n* duj 



i=T,L 



E d. 

i=T,L 

327r2 
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2 2uo 
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2 , Ji' 

7 + + 
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27r 
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1 + 2n,)^i - B 
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In 2 



ln^2 



(51) 



with ipT^L = Disc(A^_^nf,_^) and 
E diBr' = -d 



IT 



'Mi 
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Ml 
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